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PREFACE 


The High-Speed Research Program sponsored the NASA High-Speed Research Program 
Aerodynamic Performance Review on February 8-12, 1999 in Anaheim, California. The review 
was designed to bring together NASA and industry High-Speed Civil Transport (HSCT) 
Aerodynamic Performance technology development participants in areas of: Configuration 
Aerodynamics (transonic and supersonic cruise drag prediction and minimization) and High-Lift. 
The review objectives were to: (1) report the progress and status of HSCT aerodynamic 
performance technology development; (2) disseminate this technology within the appropriate 
technical communities; and (3) promote synergy among the scientist and engineers working HSCT 
aerodynamics. The HSR AP Technical Review was held simultaneously with the annual review of 
the following airframe technology areas: Materials and Structures, Environmental Impact, Flight 
Deck, and Technology Integration. Thus, a fourth objective of the Review was to promote synergy 
between the Aerodynamic Performance technology area and the other technology areas within the 
airframe element of the HSR Program. 


The work performed in the Configuration Aerodynamics (CA) element of the High-Speed 
Research Program during 1998 was presented in the following sessions: 


Propulsion Integration 
Analysis Methods 
Design Optimization 
Testing 


The work performed in the High Lift (HL) element of the High-Speed Research Program during 
1998 was presented in the following sessions: 


High-Lift Configuration Development 
Tools and Methods Development 


The proceedings for the Aerodynamic Performance Annual Review are published in three 
volumes: 


Volume I, Parts 1 and 2 Configuration Aerodynamics 
Volume II, Parts 1 and 2 High Lift 


AP Review Chairperson: David Hahne 
NASA Langley Research Center 
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LaRC conducted a code validation study for the OVERFLOW code to ascertain its accuracy 
for boattail drag prediction. The OVERFLOW results compared favorably with the LaRC 
16-ft TWT data, and pnor CFD solutions from PAB3D and CFL3D. 


The ultimate goal ts to investigate the installation drag of the nacelle boattails with powered 
nozzles at transonic Mach numbers. The OVERFLOW solver was chosen because of its 
ability to accept volume overlapping structured grid for very complex airframe 
configurations. Structured grid components for representing the transonic nozzle boattail can 
be added to the BCAG grid for a TCA airframe with 2D bifurcated inlet and flow through 
nacelle without alteration.. The focus of this research was to determine the suitability of the 
OVERFLOW solver for accomplishing this ultimate goal. 


Configuration Aerodynamics Technology Development 
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This study was performed under High Speed Research II - Configuration Aerodynamics 
WBS 4.3.1.3.2: Propulsion Induced Effects 

Milestone 4.3: Initial Propulsion Induced Effects 

Coordination Responsibility: NASA Ames Research Center 
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@ N1-S1 Rectangular Nozzle Boattail Geometry 


@ Review of Aerodynamic Features and Comparisons to 
Prior PAB3D Solutions 


@ Comparison of OVERFLOW Solutions to Experiment and 


Other CFD Solutions at M= 0.9, 1.11, 1.20 
@ Boattail Drag Prediction Assessment 
@ Work Plan for FY99 





This presentation will first introduce the transonic nozzle boattail wind-tunnel model 
geometry, followed by an examination of aerodynamic features based on the current 
OVERFLOW solutions and the solutions obtained previously using PAB3D, comparisons of 
Cp on the flap surface between the OVERFLOW solutions, wind tunnel data, and solutions 
from other CFD codes, an assessment of boattail drag count prediction, and a work plan for 
FY99. 


Isolated Nozzle Geometry (Transonic) 
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The N1 flap has a circular are profile between its hinge line and the flap trailmg edge at the 
nozzle exit. The nominal nozzle flap angle setting is 8 degrees in reference to the nozzle 
horizontal plane of symmetry. Surface slope at the trailing edge is actually 16 degrees. The 
S1 sidewall has a taper angle of 4 degrees. The top of the side wall is rounded by a quarter 
circle according to the local thickness. The computational mesh represents a quarter of the 
nozzle configuration with two planes of symmetry. The upstream portion of this isolated 
nacelle model, not show in this figure, is a pointed forebody with super-elliptic cross 
sections. The total model length is 64 inches. The computational mesh used here for 
OVERFLOW analysis, other than volume overlap between blocks, is identical to the 
multiblock grid used for PAB3D (LaRC) and CFL3D (Boeing-LB) computations. The 
nozzle internal convergent-divergent flow path and the supersonic jet exhaust flow following 
the nozzle exit are included in the computational domain. 


Mach Number Contours at Mach = 0.9 
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Two versions of the OVERFLOW solutions are shown on the left hand side: one based on 
the SA turbulence model and another based on the Menter’s SST model. A corresponding 
PAB3D solution is shown on the right hand side which is computed using a two equation k-e 
turbulence model and the Girimaji algebraic Reynolds stresses. A local transonic expansion- 
compression pressure pattern can be seen in all three solutions. A small region of flow 
separation can be seen near the flap trailing edge. The nozzle internal flow path and the 
exhaust flow are also shown in this figure. The nozzle pressure ratio (NPR) is 5.0 for all 
three cases: M=0.9, 1.11, and 1.20. Important flow features to be demonstrated in the next 
three sets of flow quality figures are: flow expansion and recompression over the nozzle flap, 
shock boundary layer interaction, flow separation, and the formation of vortices as a result of 
cross flow over the top of the side wall. 


Normalized Pressure Contours on Nozzle Flap 
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An isometric view of the normalized contour over the nozzle flap is shown mm this figure. 
The pressure pattern near the hinge line is three dimensional m nature. As the flow proceed 
down the flap, the confinement of the side wall causes the pressure distribution to show 
typical two dimensional features. A particle trace, in red color, shows the position and width 
of the induced vortex over the side wall. 


Mach Number Contours at M = 0.9 
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Mach number contours for a cross section at the trailing edge of the flap are show for the 
same solutions at M=0.9. The flow outside of the side wall is pulled over the top of the side 
wall toward the low pressure region above the nozzle flap. A strong vortex is formed 
inboard of the side wail. The induced change in pressure distribution on the flap surface 
should reflect a drag component equal to the momentum loss carried away by the detached 
vortex element. The difference in thickness of flow separation at the trailing edge of the flap 
can be seen in this figure. 


Mach Number Contours at M = 1.11 
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The OVERFLOW solutions for M=1.1]1 was computed by using the SA-turbulence model 
alone. The flow expansion over the flap starts at the hinge line, and continues much farther 
downstream over the flap. The expansion was eventually termmated by a shock. Shock 
boundary layer interaction has crated a flow separation which extend from the root of the 
shock down to the flap trailing edge. Since the NPR for the internal flow remains the same 
as the M=0.9 case, the internal and exhaust flow characteristics remains the same as before. 


Normalized Pressure Contours on Nozzle Flap 
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The normalized pressure contours for the OVERFLOW and the PAB3D solutions show 
significant differences in pressure distribution of the flow surface where most of the drag 
force is produced. While the low pressure area of the OVERFLOW solution covers a much 
larger area, the pressurization near the trailing edge is also much stronger. The particle 
traces for the two solutions are also quite different. The vortex as shown by the 
OVERFLOW solution is more diffused in appearance. Furthermore, clear entramment paths 
from a significant length of the side wall directly to the vortex core region. The PAB3D 
solution, on the other hand, does not show such a pattern. 


Mach Number Contours at M = 1.11 
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Cross Section View at X = 64.0 inch 
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The vortex strength as shown by the OVERFLOW solution appears to be much more 
diffused than the corresponding vortex structure in the PAB3D solution. Flow separation 
over the flap surface for the OVERFLOW solution appears to be much thinner than its 
counterpart in the PAB3D solution. However, these qualitative flow features do not tell 
much about the change m drag count. 
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Velocity Vector Distribution at M = 1.11 
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Cross Section View at X = 64.0 inch 
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OVERFLOW: SA Turbulence Model PAB3D: ASM Turbulence Model 


In order to provide a better understanding of the quantitative nature of the vortices depicted 
in the OVERFLOW and the PAB3D solutions, velocity vector plots are shown in this figure. 
Velocity vectors of the same magnitude are located farther away from the vortex core inthe 
OVERFLOW solution versus the PAB3D solution. Again, the significance of these patterns 
in terms of drag forces is not obvious with a quantitative analysis of various aerodynamic 
quantities such as mass flux over the side wall, momentum flux, and total circulation around 
the vortex. 
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Mach Number Contours at M = 1.2 
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OVERFLOW: SA Turbulence Model PAB3D: ASM Turbulence Model 





The OVERFLOW and the PAB3D solution are very similar inthis comparison. There are 
differences in shock position as well as the size and thickness of flow separation near the 
trailing edge of the nozzle flap. The shock is definitely stronger than the M=1.11 case, and 
the shock position is further downstraam on the flap surface. Again, the internal and the 
exhaust flow regions remain similar to the M=0.9 and 1.11 cases since the NPR remains the 
same as before: NPR = 5.0. 


12 


Normalized Pressure Contours on Nozzle Flap 
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Other than the strong influences of the vortex, the flow expansion and compression patterns 
on the flap surface is decidedly two dimensional in appearance. As was in the M=1.11 case, 
the OVERFLOW solution, in comparison to the PAB3D solution as shown on the right hand 
side, has larger regions of low pressure as well as stronger recompression hear the trailing 
edge. In essence, these two effects cancels each other as far as overall drag forces are 
concern. 
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Mach Number Contours at M = 1.20 
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OVERFLOW: SA Turbulence Model PAB3D: ASM Turbulence Model 





Similar to the M=1.11 case, the PAB3D vortex is much stronger than the vortex in the 
OVERFLOW solution. The influence of the vortex on the flow field above the flap and near 
the side wall is evident in this figure. 


Comparison of Surface Pressure Coefficients 


" AE RRS AE WRIT NTE, G2 php Sty Sr ahe 


Center Line: y = 0. Quarter Span: y = 1.31 


EFO Top Cantetiine yod in 


‘a+ QHERFLOW, Monte $37 un. CBLIO, SpetatAfimerse 
-- OFLIO, Spatert-Alnersa 


EFD Top Outhaoré yu2.¢1 in 
'- OVERFLOW, Spatart-Almares 
--~ QVERFLOW, Menter SST 
CFL3O, Spatert-Alename 





Wind tunnel measurement data and CFD results from OVERFLOW, PAB3D, and CFL3D 
are shown in this figure. The CFL3D solutions were obtained by Boeing Phantom Works at 
Long Beach. At M=0.9, the flow initiated supercritical expansion near the hinge line, and 
continued over the initial portion of the flap surface. Pressure recovery started at 
approximately at model station 57. Two OVERFLOW solutions, using the SA and the SST 
turbulence models are shown m this figure. 


Good agreement is shown comparmg to wind tunnel data and pressure coefficients computed 
by other codes. The OVERFLOW solutions show a stronger and more extended flow 
expansion over a significant portion of the nozzle flap. The pressure recovery of the 
OVERFLOW solution achieved a higher pressure at the trailing edge of the flap in 
comparison to the data and the PAB3D and CFL3D solutions. 
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Comparison of Surface Pressure Coefficients 
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The nozzle boattail measurements conducted in the LaRC 16-FT transonic wind tunnel had 
produced an accurate database for boattail drag. However, the overall model length of 64 
inches was long than most other propulsion component models tested in this tunnel. At 
M=1.11, shock reflected from the tunnel wall had pressurized the base of the nozzle boattail. 
As a result, the measured pressure coefficient was raised to much higher values, with the 
strongest influence near the trailing edge of the nozzle flap. The correct value for pressure 
coefficient at this mach number is much better represented by the PAB3D solution. 


Hence, only the CFD solutions between OVERFLOW and PAB3D can be compared. Flow 
expansion between the hinge line and model station 61 as predicted by the two solutions are 
nearly identical. Flow expansion in the PAB3D solution is termmated by a shock shortly 
downstream of station 61. The OVERFLOW solution shows further expansion before a 
shock is set up farther downstream. As was in the M=0.9 case, the recovery pressure at the 
flap trailing edge is higher for the OVERFLOW solution. 
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Comparison of Surface Pressure Coefficients 
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Wind tunnel measurement data and CFD results from OVERFLOW, PAB3D, and CFL3D 
are again shown in this figure. The CFL3D solutions were obtained by Boeing Phantom 
Works at Long Beach. At M=1.2, the pressure coefficients throughout the expansion phase 
up to model station 61. Although both OVERFLOW and CFL3D solutions were computed 
by using the SA turbulence model, the shock positions are downstream and upstream of the 
data for the two solutions, respectively. Both OVERFLOW and CFL3D solutions show 
recovery pressures which are higher than the measured data near the flap trailing edge. From 
previous studies of nozzle boattail pressure recovery using the PAB3D code, turbulence 
models have a strong effect on pressure recovery characteristics in the transonic speed range. 
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Integrated Boattail Drag Summary 
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Cdp Cd 
Experiment PAB3D OVERFLOW Experiment PABID OVERFLOW 


0.000065 0.000057 0.000502 SA 0.000137 0.000507 
0.000503 SST 0.000508 
Flap 0.000109 0,000117 0.001207 SA 0.000187 0.001210 
0.001211 SST 0.001214 
Nozzle 0.000174 0.000174 0.001710 SA 0.000324 0.001717 
0.001715 SST 0.001721 
M=1.11 
Sidewall 0.000062 0.000295 SA 0.000257 0.000299 
Flap 0.000289 0.000589 SA 0.000864 0.000593 
Nozzle 0.000351 0.000884 SA 0.001121 0.000892 
M=1.20 
Sidewal! 0.000168 0.000253 SA 0.000211 0.000258 
Flap 0.000843 0.000473 SA 0.000882 0.000467 
Nozzle 0.001011 0.000726 SA 0.001161 0.001093 0.000734 





These numbers shows a complex picture for boattail drag predictions in the transonic range. 
According the IMS method proposed by Hoyt Wallace (BPW), the presence of side walls 
would significantly modify the Mach number trends of nozzle boattail drag through the 
transonic range. The drag rise begins at approximately M=0.9. Drag rises rapidly bey ond 
=0.95, and contmue to rise through M=1.1. For nozzle boattails without side walls, the 
drag coefficient begins to decline rapidly beyond M=1.10. However, the drag coefficient 
would maintain its peak value through M=1.2 before declining to much lower values. 


The nozzle boattail measurements conducted in the LaRC 16-FT transonic wind tunnel! had 
produced an accurate database for boattail drag. However, the overall model length of 64 
inches was long than most other propulsion component models tested in this tunnel. At 
M=1.11, shock reflected from the tunnel wall had pressurized the base of the nozzle boattail. 
As a result, the measured drag coefficient was less than half of what would be expected by 
using the IMS method. The correct value for drag coefficient at this mach number is much 
better represented by the PAB3D solution. 


Drag count computed using FOMOCO and the OVERFLOW solution at M=0.9 is 17.2, 
much higher than anticipated from either the experiment or the PAB3D solution. Drag 
counts at M=1.11 and 1.20 are 8.92 and 7.34, respectively , somewhat lower than the 
measured values for this configuration. The source of the differences are being investigated. 
The most likely sources of variation are the difference between upwind and central 
difference schemes and differences in turbulence models used in these computations. 

How ever, it remains possible that procedural errors could have been made in the force 
integration process. 

Most of the drag forces resulted from flow expansion over the deflected nozzle flap. The 
side wall drag contribution is small for this geometry because the taper is 4 degrees. The full 
scale nozzle design would require a taper of at least 8 degree for closure of a much thicker 
side wall. Hence, the differences in side wall drag count between awindt. 1 model and 
the full scale design should be noted. 
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Overall Remarks 





“rn 


EES 


@ OVERFLOW VALIDATION 
Flow feature predictions are similar to prior PAB3D solutions 
» Cp distribution compared favorably with LaRC 16-FT wind tunnel 
data and PAB3D solutions 
Shock position and magnitude of pressure recovery are different 
between CFD solutions and measured data (CFL3D solution was 
obtained by BWP-LB) 
Variation of drag count predictions is noted in this study; 
numerical scheme, turbulence model, or simply procedural error in 
force integration may be responsible 
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@ Aerodynamics 

» Understanding of aerodynamic features can lead to nozzle boattail 
design improvements and drag reduction 

» CFD drag prediction is a useful tool for expanding nozzle boattail 

drag data base beyond range of wind tunnel measurements 








No notes for this slide. 


Future Plans 
e HSR - FY99 a 


» Add transonic nacelle grid components to TCA chimera grid 
constructed by BCAG (Chaney) 


» Compute installed nozzle boattail drag components using 
OVERFLOW with SA turbulence model 


» Powered effect will be included in the analysis 


@ Computational Methods Development 
» Add chimera grid capability to PAB3D 
» Compute same cases and compare all available results 





The installed nacelle boattail study using OVERFLOW and amodified TCA chimera grid for 
Mach numbers 0.9, 1.10 and 1.20 will be completed within FY99. However, support for 
further studies of the installation effect, such as change in geometry or drag optimization will 
not be available. The proposal to extend PAB3D for chimera capability is based on its 
excellent ability to predict propulsion component performance for both intemal! flow and 
configuration aerodynamics. Support for such work could come from the base program or 
other sources. 
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High Speed Configuration Aerodynamics 


Isolated PIE Nacelle 





¢ 2-D Inlet and Ramp, 2-D Nozzle (aspect ratio 1.5) 
— More complicated in CFD grid generation than Ref H 2-D 


e Operational and Aero Reference Nozzles 
¢ Mach 0.9 and 1.1 Nozzle Configurations 


¢ Nozzle Flaps Longer than Side-walls 
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Analytic Methods Test Programs 
Approacnes and Applications and Techniques 


Methods Down Select Nacelle / Diverter WT Database 
—<——$—$__—_——] | | Design Integration 
Program “Viscous Drag Prediction © WT Data Corrections 
; ————s Tech. Baseline : ; 
Cruise Point Optimization Development High Re. No. Testing 
Multi-Point Optimization Aero S&C PIE Test Program 


S&C CFD Predictions Development 





The assessment of the CFD flat plate fully turbulent flow skin friction 
predictions is an element of the “Viscous Drag Prediction” technology 
development element shown in the Configuration Aerodynamics 
program on a page shown in this figure. 
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Topics 


Variations in Viscous Drag Predictions 
Resolution Strategy 

Flat Plate Skin Friction Comparison Database 
BLB-PW CFL3D Predictions 

BCA Overflow Predictions 

Ames OVERFLOW Predictions 


Conclusions 





Recent CFD validation studies have shown significant variations in viscous drag predictions 
between the various methods used by the NASA and industry HSCT organizations. The 
methods include Navier Stokes CFD codes in which the viscous forces are part of the 
solutions, and predictions obtained from the different fully turbulent flow flat plate skin 
friction drag equations used by the various organizations. 


In this paper, the variation of these viscous drag predictions will be shown. The strategy 
developed to resolve these differences will be discussed. The first step in the resolution 
Strategy was the development of a skin friction database for flat plate fully turbulent flow. 
This database will be briefly reviewed. The comparisons of CFD skin friction prediction by 
Boeing Phantom Works Long Beach - , BPW-LB, by Boeing Commercial Aircraft 
aerodynamics, BCA, and by NASA Ames will be reviewed. 


The study conclusions will be summarized. 
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Comparison of Predicted Viscous Drags From Various 
CFD Codes for the TCA Wing/Body Configuration 


M.. = 2.4 ,C, =0.1, Re = 6.36x10° 


CFD Analyses 
AIRPLANE - Nicolai 
2 SYNB7MS - Neolai 
3 GFLSD - Van Onest 
TRANAIR - Sommer-Short 


Ormax 

5 OVERFLOW - Spalat-Alimaras- Ames 
OVERFLOW - Spalart-Allmaras ~ BCA 
UPS - Baldwin-4.omax 


¥ 
S 
= 
o 
rx) 
> 
=) 
© 


CFD Code 





This figure illustrates the variations in viscous drag predictions for the 
TCA wind tunnel model wing plus body configuration. There are 
significant variations in flat plate theory predictions used in the inviscid 


CFD analyses as well as the CFD predictions obtained with the viscous 
analyses. 
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Comparison of Predicted Viscous Drags From Various CED 
Codes for the TCA Wing/Body/Nacelle/Diverter Configuration 


M..=2.4,C, = 0.1, Re =6.36x108 





This figure shows similiar comparisons for the wing plus body plus 
nacelle configuration. There is seen to be a three drag count variation in 
the flat plate predictions corresponding to the inviscid analyses, 
(numbers 1 through 4). The BLB-PW CFD predictions using CFL3D 
appear to match the flat plate theory predictions. The OVERFLOW 
predictions of BCA and Ames predict significantly lower drag levels. 
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Resolution Strategy 


ae Bp * Establish Data for Validation of Skin Friction Predictions 
-. Database. and Predictions * Evaluate Flat Plate Skin Friction Theories 


sessment of CFD Predictions * Evaluate CFD Predictions of Fully Turbulent Skin Friction 
of Flat Plate Skin Friction Select Turbulence Model (s) for subsequent analyses 


Symmetric Mode! Wind tunnel => ¢ Obtain Drag Data at Zero Lift 
Test Programs ¢ Realistic Configuration Geometry 


sap > Database for Validation of Prediction of Mach Effects 
“hitestRogam | same + Database for Validation of Prediction of Rel Effects 


Assessment of CFD Predictions 
of Symmetric Mode! 


emp + Validate CFD Predictions of Mach Effects 
[NTP tesrPrecietens Analyses | maa « Validate CFD Predictions of Reynolds Number Effects 


mee + Evaluate Refine CFD Calculation Processes 


Biol CFD Precictions ¢ Determine if Closure Has Been Achieved in Viscous 
pind Drag Predictions 





This illustrates the strategy that was developed to resolve the viscous drag 
prediction differences. This consisted of a series of sequential activities: 


° Establish a database of fully turbulent flow flat plate skin friction data to be 
used for the validation of the corresponding CFD predictions. Flat plate 
prediction methods were assessed. A modified flat plate skin friction prediction 
method was developed that accurately represents the mean of the test data and 
captures both the Reynolds number and Mach number variations of this mean. 


¢ The second step includes the comparison of the CFD predictions of fully 
turbulent flow flat plat viscous drag with the mofified flat plate theory. The results 
of this activity is the subject of this paper. 


¢ A symmetric model representation of the HSR TCA configuration was defined 
and will be fabricated and tested to obtain data for validation of CFD viscous 
drag predictions on an HSR type configuration, Supersonic tests are planned in 
the Boeing Polysonic Wind Tunnel to obtain supersonic drag data at moderate 
Reynolds numbers The model was also planned to be tested in the NASA . 
Langley NTF tunnel to obtain data for a wide range of Reynolds numbers. 


¢ The final element is to recalculate the drag of the TCA to see if the variations 
between the theories has vanished and the theory predicts the test results. 


eThe elements that are crossed out have been canceled by reduction in program 
funding 
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Why the Interest in Flat Plate Turbulent Boundary layers ? 






¢ First Step in Evaluating Navier Stokes Prediction Methods 






¢ Help Sort Out Appropriate Turbulence Models 





* Good Estimate of Viscous Drag of HSCT Type Configurations 
( Easy, Quick , Robust and Accurate) 













¢ PD Drag Prediction Methods 


¢ Extrapolation of Wind Tunnel Data to Flight Conditions 





* 6 Predictions Used to Size Diverter Height 






¢ 6* plus CF Predictions Used to Calculate Spillage and Internal 
Drag of Flow-Through Nacelles 





* Quick Estimate of Surface Temperature 






¢ Provides Physical Insight into Viscous Flow Characteristics 






It is felt that the first step in validating the viscous drag predictions of any Navier Stokes 
code is to make sure that predictions of the local and average skin friction drag and 
boundary layer must match the “simple” flat plate measured test data over the range of 
Mach numbers and Reynolds for which the codes will be used. This process will help to 
evaluate the applicability of of the various turbulence models. 


Because HSCT configurations have rather thin wings, slender bodies and low cruise lift 
coefficients, experience has shown that flat plate skin friction calculations provide good 
estimates of the viscous drag of HSCT type configurations. The predictions are easy, quick, 
robust and quite accurate. 


The current PD viscous drag prediction methods are based on flat plate skin friction drag 
calculations. Currently wind tunnel data is extrapolated to flight conditions using flat plate 
friction drag predictions. 

Flat plate estimates of the boundary layer thickness are used as the preliminary criteria for 
specifying the boundary layer diverter height for the HSCT nacelle installations. Boundary 
layer displacement thickness predictions together with CF calculations are used to calculate 
the spillage and internal drag of wind tunnel flow though nacelles. 

Local skin friction calculations corrected for local dynamic pressure effects can be used to 
estimate local surface temperatures. 


The boundary layer thickness information presented in this note also provides some 
physical insight in to the fundamental features of turbulent flat plate flow. 
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Incompressible Local Skin Friction Data 
[ ACT = Crest - Chrneory ] x 104 [ ACT / Chingory ] x 10? 


« Theory: Modified Schultz / Grunow 
0005 |° Test Data: Kulfan, R. M.. D6-7161 
Rex 108 


beret caer se 3 are +16 
! ae +o, PFAe- Mean 
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Flat plate skin friction data was obtained from a number of experimental sources. These 
data cover a wide range of Mach numbers and Reynolds numbers. Comparisons were 
made with various flat plate theories to select the theory that most closely matches the test 
data. The results of these assessments are presented in the Reference shown below. 


The flat plate theories are based on the reference temperature method. This method 
assumes that the incompressible skin friction equations apply to supersonic Mach numbers 
provided that the density and viscosity are calaculate at some reference temperature that 
represents the variation of temperature across the boundary layer. 


This figure shows the comparison of the modified Shultz / Grunow equation with 
incompressible test data. Statistical analysis of the differences between the test data and 
corresponding Cf predictions shows that the mean of the differences is ACf = -.000000671 
which corresponds to an average difference of 0.13% .The standard deviation of data about 
the mean is approximately 0.7 counts of drag ( ACf = 0.000067) which corresponds to 2.8% 
of the corresponding predicted value. 


The modified Shultz / Grunow equation therefore appears to provide an accurate estimate 
of incompressible local skin friction coefficient over the entire range of Reynolds Numbers 
covered by the test data. 


Reference: Kulfan, R. M., “Historic background on fiat plate turbulent flow skin friction and boundary layer growth”, HSR Airframe 
Technical Review, Feb 1998 
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Comparison of Compressibility Effects Predictions 
Skin Friction Ratio Cf/Cfi CF/CFi 








Skin Friction Ratio Cf/Cfi CF/CFi 






Mach Number 






=~ Theory: Monaghan Rex = 5x10 

——— Theory: Sommer / Short Rex = 5x106 

© Test Data: Kulfan, R. M. D6-7161 

© Test Data: White, F. M. , “Viscous Fluid Flow" 






This figures shows some of the compressible flow skin friction datya used to validate the 
flatplate theories. This compares the compressible skin friction predictions obtained using 
two commonly used T* methods, the Monaghan T* and the Sommer-Short T* method. 


The Sommer-Short T* €quation results in compressible skin friction values consistently 
higher than predicted using the Monaghan method. It was for this reason that the Boeing 
US SST program switched from the Monaghan method to the Sommer-Short method. 


The full scale SST performance predictions were obtained from wind tunnel data corrected 
to full scale conditions. Wind tunnel skin friction drag is higher than the full scale conditions. 
Using higher skin friction values calculated by the Sommer -Short method resulted larger 
Skin friction corrections. This resulted in higher L/D assessments for the SST. 
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Evaluation of Reference Temperature Equations 
Monagham Sommer and Short Kulfan 
T/T, =1 + 0.1246 M,? T/T, =1 +.0.1151 M,? T/T, =1 +0.1198 M,? 
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Statistical analyses were made of the differences between Cf predictions and the corresponding test 
data as shown in this figure . The theoretical predictions were obtained using three different T* 
equations. 

The “‘scatter” in the test - theory increments are essentially equal. The mean of the differences 
between the test and theory, however differs between the predictions obtained using the different T* 
equations. 


The “mean” of the theory - test differences obtained using the Monaghan T* equation is 
approximately 1% low. The “mean” of the theory - test differences obtained using the Sommer- 
Short T* equation is approximately 1% high. The constant for the Kulfan T* equation was therefore 
chosen to be the average of the Sommer-Short and the Monaghan constants. 


This essentially resulted in a mean error between the test data and the theoretical predictions of zero. 


The test data scatter about the mean has a standard deviation of about 4.5%. This large scatter is in 
part due to the variations of Reynolds number of the test data. The Reynolds number for the test data 
10° to 10’. 
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Conversion of Compressible Cf Data to E 
* Kulfan T* Method 
¢ Modified Schultz-Grunow Cf Equation 


quivalent Incompressible (Cfi)., Data 






Cf 





* (Rex), z (Rex) x (T/T* pfu") re 
° (Cfi) og = (T*/T) xCf 






© Mach 0.20 Matting, F. W., et al 
© Mach 2.00 Shutts, Ww. H., etal a as 
~ X Mach 2.95 Matting, F. W., et al mein 
| : | Re pas os! 
§ 7 7 
10 Rex 10 108 (Rex), 10 108 


Reynolds Number Equivalent Incompressible Reynolds Number 















The T* equations can also be used to convert the compressible skin friction to equalivent 
incompressible data. This transformation procedure, as shown in the Figure, “collapses” all 
of the test data about the incompressible skin friction curve. This approach can provide a 
convenient means to assess the accuracy of the theoretical methods to account for 
compressibility effects simultaneously over a range of Mach numbers and Reynolds 
numbers. 
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Conclusions - Flat Plate Database Developemnt 


© Modified Incompressible Equations and Improved T*/T Method 
Predict “Mean” of Available Flat Plate Skin Friction Drag Measurements 


e New Methods Presented That Appear to Provide Good Estimates 
of Boundary layer Thickness and Displacement Thickness 


¢ Compressibility Effects Have Very Little Effect on The Shape or Height 
of the Turbulent Flat Plate Velocity Profile. 


e Boundary Layer Displacement Thickness Increases Rapidly With Mach Number 


° Comparisons of Navier Stokes CFD Predictions of Flat Plate Turbulent Skin 
Friction Drag and Boundary Layer Growth, With the Test Data and / or Theory 
Presented in This paper, is considered to be a Necessary and Vital Step to 
Validating the Codes For HSCT Viscous Drag Predictions . 


* Need Additional / Quality Experimental CF Data. 
- Locate Available Existing Data 
- Symmetric Model Tests 
- Segmented Axisymmetry Body of Revolution 
- Utilize TU-144 Flight Test Data 
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The modified incompressible CF equations and the improved T* equation presented in the 
reference paper appeared to consistently match the test better than the other flat plate CF 
methods currently in use on the HSCT program. It was recommended that the methods 
presented there, be adapted as the official HSCT flat plate calculation methods. 

The boundary layer thickness, and displacement thickness calculations methods presented 
in that paper seem to be validated by the existing data. 

Compressibility effects were shown to have little effect on either the shape or height of a 
turbulent boundary layer. The displacement thickness however varies rapidly with 
increasing Mach number. | 

A modified Shultz / Grunow incompressible local skin friction equation and the modified 
Prandtl / Schlichting average skin friction were used with the Kulfan T* equation in the 
studies reported in this paper to evaluate the CFD predictions of fully turbulent flow flat plate 
skin friction drag. 
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BPW-LB Skin Friction Analyses Analyses 





* Code: CFL3D 







Average Skin Friction, CF 





¢ Turbulence Models: 
* Baldwin - Lomax 

* Spalart - Allmaras 
“ Menters SST 









Mach Numbers: 
"OS 
"15 
"2.25 
"25 










Reynolds Number: 
* 10° to 200 x 108 





The BPW-LB average skin friction predictions were made using CFL3D 
and a number of turbulence models fora range of Mach numbers and 
Reynolds as shown in the figure. 
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CED Calculated Average CF Comparisons with Flat Plate CF 
« BPW-LB CFL3D Code 
¢ Baldwin - Lomax Turbulence Model 








BERBERBSERESES 


Average Skin Friction , CF x 10* 


~~ — at = = FA) N f%& 
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107 
Reynoids Number, Rex 


This compares average skin friction predictions obtained using the 
Baldwin - Lomax turbulence model, with the flat plate predictions. The 
Calculations were made for Mach = 0.5, 1.5, 2.25 and 2.5 
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Differences Between CFD Calculated CF and Flat Plate CF 
¢ BPW-LB CFL3D Code 
* Baldwin - Lomax Turbulence Model 


Differences in CF Predictions: CF een -CF ep Differences in CF Predictions: (CF cep -CF pp V/CFep 


-- Percent 


(CFegp -CF ep) x 1 Or4 


(CF cep -CFepV/CF pp 
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10” 108 10? 108 
Reynoids Number, Re! Reynolds Number, Rel 
CFD Calculated CF* for Mach e 0.50 
CFD Calculated CF* for Mach = 1.50 


CFO Calculated CF* for Mach = 2.25 
CFD Calculated CF* for Mach = 2.50 


3 





This shows the differences between the CFL3D predictions and the flat 


plate theory both as incremental differences, and differences in percent. 


The differences in the predictions are quite Mach number dependent. 
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Equivalent incompressible Average Skin Friction Calculations 
¢ BPW-LB CFL3D Code 
¢ Baldwin - Lomax Turbulence Mode! 


50 
48 
46 | 
44 
42 
40 |" 
38 
36 
3a |} 
32 
30 | 
23 |: 
26 
24 
22 
20 | 


2 
<x 
= 
2 
a 
= 
s 
8 
L. 
O 


for Mach = 0.00 
CFO Calculated CF* for Mach = 0.50 
CFD Calculated CF for Mach = 1.50 
CFD Calculated CF” for Mach = 2.25 
CFD Calculated CF* for Mach = 2.50 
Mean of CFD Calculations of CF" 





10’ 
Equivalent incompressible Reynolds Number, Rex 





The Kulfan T* equation was used to transform the CFL3D predictions to 
incompressible skin friction data. The dash red line in this picture is the 
mean of the CFL3D predictions. It appears that the CFL3D predictions 
with the Baldwin-Lomax turbulence matched the flat plate predictions 
and the variation of Cf with Reynolds number 


268 


Average Skin Friction Calculations: CF /CFi 


« BPW-LB CFL3D Code 
° Baldwin - Lomax Turbulence Model 


CF/CFi 


Fiat Plate CF for Rel = 10°6 
‘| Flat Plate CF for Rel = 2x108 
CFD cakulated CF for Rai = 10°G 
CFD calculated CF for Rel = 2x10“8 
CFD calculated CF for Ref = 106 to 2x10°8 


0.2 04 0.6 0.8 1.0 4.2 1.4 
Mach Number 
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CFD Calculated Average CF Comparisons with Flat Plate CF 
- BPW-LB CFL3D Code 
e Spalart - Allmaras Turbulence Model 





Average Skin Friction , CF x 104 


108 10’ 10° 


Reynoids Number, Rex 
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Differences Between CFD Calculated CF and Flat Plate CF 
¢ BPW-LB CFL3D Code 
* Spalart - Allmaras Turbulence Mode! 


Differences in CF Predictions: CFegp -CF ep Differences in CF Predictions: (CFegp -CFep/CF ep 


Percent 


(CF erp -CF ep)x 1044 
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CFD Calculated CF” for Mach = 0.50 
CFO Calculated CF* for Mach = 1.50 
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CFD Calculated CF* for Mach = 2.50 
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Equivalent incompressible Average Skin Friction Calculations 
¢ BPW-LB CFL3D Code 
¢ Spalart - Allmaras Turbulence Model 
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CFD Caiculated CF” for Mach = 0.50 
CFD Calcutated CF* for Mach = 1.50 
CFD Calcutated CF* for Mach = 2.25 
CFD Calculated CF* tor Mach = 2.50 
Mean of CFD Calculations of CF° 
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Average Skin Friction Calculations: CF /CFi 
* BPW-LB CFL3D Code 
* Spalart - Alimaras Turbulence Mode! 







Bes ME aie Renenne SHENAE es ahs 4 Sees Gaede doi ssecavaas vslneteneeineaiic’es Grdh Laced amishescon cp isl 


CF/CFi 







Flat Plate CF for Rel = 10°6 
‘| Flat Plate CF for Re! = 2x17%8 

CFD calculated CF for Rel = 10°6 

CFD caicutated CF for Re! = 2x10°8 

CFO calculated CF for Rel = 1046 jo 2x18 
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CFD Calculated Average CF Comparisons with Flat Plate CF 
¢ BPW-LB CFL3D Code 
¢ Menters Turbulence Model 


Average Skin Friction , CF x 10* 
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Differences Between CFD Calculated CF and Flat Plate CF 
¢ BPW-LB CFL3D Code 
¢ Menters Turbulence Model 


Differences in CF Predictions: CFery -CF rp Differences in CF Predictions: (CFegy -CFrpV/CF ep 
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Equivalent incompressible Average Skin Friction Calculations 
¢ BPW-LB CFL3D Code 
¢ Menters Turbulence Model 
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Average Skin Friction Calculations: CF /CFi 


¢ BPW-LB CFL3D Code 
* Menters Turbulence Modei 


Flat Plate CF tor Rel < 10°6 
‘| Flat Plate CF for Rel = 2x10°8 
CFO cakulated CF for Rel = 106 
‘| CFD cakulated CF tor Ael = 2x10°9 
CFD calculated CF for Rel = 106 to 2x18 
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BLB CFL3D Flat Plate Viscous Drag vs Fiat Plate Skin Friction Drag Calculations 


Ditterences in Percent 
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BAC Skin Friction Analyses Analyses 


¢ Code: OVERFLOW 
¢ Local Skin Friction, Cf 
* Average Skin Friction, CF 


e Turbulence Models: 
* Baldwin - Lomax 
* Spalart - Alimaras 
*Menters SST 
* Baldwin - Barth 
*k-e 
*"kK-Ww 


e Grid vertical spacing variations 


¢ Mach Numbers: 
"0.9 
"2.4 


¢ Reynolds Number: 
* 105 to 6 x 106 
* 105 to 200 x 10° 
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BCA OVERFLOW Local Skin Friction Calculations, Cf 


Mach = 0.9: Baldwin - Lomax Turbulence Model 
Unform vertical grid spacing 
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BCA OVERFLOW Local Skin Friction Calculations, Cf 


Mach = 2.4: Bladwin - Lomax Turbulence Model 
Unform vertical grid spacing 
: 2 ee 
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BCA OVERFLOW Local Skin friction Calculations Comparisons with Flat Plate Cf 
Baldwin - Lomax Turbulence Model 

Unform vertical grid spacing 
Local Skin Friction Calculations 2 ee 
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Difference in Cf predictions - Percent 
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BAC OVERFLOW Average Skin friction Caiculations Comparisons with Flat Plate CF 
Baidwin - Lomax Turbulence Mode! 
Untorm vertical grid specing 


Average Skia Friction Calculations Difference in CF predictions ~ Counts 
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BCA OVERFLOW Local Skin Friction Calculations, Cf 


Mach = 0.9: Spalart - Allmaras Turbulence Mode! 
Unform vertical grid spacing 
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BCA OVERFLOW Local Skin Friction Caiculations, Cf 


Mach = 2.4: Spalart - Allmaras Turbulence Model 
Unform vertical grid spacing 
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BCA OVERFLOW Local Skin friction Calculations Comparisons with Flat Plate Cf 
Spalart - Allmaras Turbulence Model 
Unform vertical grid spacing 
Local Skin Friction Calculations tabaci ee ee 
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BAC OVERFLOW Average Skin friction Calculations Comparisons with Fiat Plate CF 
Spalart - Allmaras Turbulence Mode! 
Untorm vertical grid spacing 


Average Skin Fricnon Calcutations Difference in CF predictions - Counts 
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BCA OVERFLOW Local Skin Friction Calculations, Cf 


Mach = 0.9: Menter's Turbulence Model 
Unform vertical grid spacing 
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BCA OVERFLOW Local Skin Friction Calculations, Cf 


Mach = 2.4: Menter’s Turbulence Mode! 
Unform vertical grid spacing 
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BCA OVERFLOW Loeai Skin friction Calculations Comparisons with Flat Plate Ct 
Menter's Turbulence Mode! 
Untorm vertical grid spacing 
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BAC OVERFLOW Average Skin friction Calculations Comparisons with Flat Plate CF 
Menter's Turbulence Model 
Unform vertical grid spacing 
Average Sidn Friction Calculations Difference in CF predictions - Counts 


St eceaeses eR ER ES 
Comme 


C8 ct - CP Sp -- 


10" 110° 


Difference in CF predictions ~ Percent 


d 
a 
‘ 
3 
if 
3 
$ 
z 
< 


110° : 110° 
Reynolds Numiscr, Rex 

— Fist Pime CF for Mach # 0.0 

=== Fist Piae CF for Mach 0.9 

— Flat Plae CF for Mach = 2.4 

aoo CFD: CF for Mach 20.9 Reb = 1.64x10%% 

eee CFD: CF for Mach = 0.9 Rel = 6.044x10% 

aes CFD. CF for Mach 22.4 Rel @ 1.97x10°6 

amex CFD: CF for Mach 2.4 Rel # 6.362 x 10%6 





292 


BCA OVERFLOW Fiat Plate Viscous Drag vs Flat Plate Skin Friction Drag 


Average Skin Friction 


Local Skin Friction 
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BCA Overflow Flat Plate Viscous Drag vs Flat Plate Skin Friction Drag 
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ARC Analyses 


Code: OVERFLOW 
Local Skin Friction, Cf 


Average Skin Friction, CF 


Turbulence Models: 
* Spalart - Allmaras 
* Menter’s SST 


Mach Numbers: 
*0.5 
"0.9 
"1.5 
* 2.0 
* 2.4 


Reynolds Number: 
* 10° to 6 x 106 
* 105 to 200 x 108 
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ARC OVERFLOW Skin Friction Calculations 
Mach = 0.5 


"Te Cf_loc: Spaiart - Allmaras 


—® Cf_loc: Menter's SST 


—#— CF_ave: Menter's SST : 
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ARC OVERFLOW Skin Friction Calculations 
Mach = 0.5 


é —— Cf_loc: Spalart - Allmaras 
—s Cf_ioc: Menter's SST : 

| ~#—CF_ave: Spatart- Alimaras [00 
~—w— CF_ave: Menter's SST : : iy pe 


1.0E+06 1.0E+07 
Rex / 1046 





297 


Mach = 2.4 


Cog Ca Baldwin - Lomax 
C221 ES) Spalart - Aimaras 
wz Ca) Menter 


Differences in percent 
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Conclusions 


* Modified incompressible CF equations and improved T*/T method 
predict “mean” of availabie fiat plate skin friction drag measurements 


¢ BPW-LB CDL3D And BCA OVERFLOW CF predictions with the Spalart - Allmaras 
Turbulence Model are consistent and show a large difference in the variation of 
CF with Reynolds Number compared to the flat plate CF predictions 


¢ The Spailart - Allmaras Cf predictions indicate a laminar to turbulent transition 
characteristic. This is the source of the large Reynolds Number variations 


* The BPW-LB CFL3D CF predictions agree better with the flatp plate theory, and 
therfore also the test data, then the BAC OVERFLOW predictions 


¢ Menter’s SST slightly underpredicts CF at low Mach Numbers, But overall 
best predicts flat plate fully turbulent flow CF 


* BAC OVERFLOW results include: 
- CFD skin friction predictions vary with vertical grid spacing scheme. 
- “constant grid spacing” near the wall, CF predictions appear better 
than “stretched grid spacing’ predictions 


* Initial assessment of ARC OVERFLOW predictions show significant differences 
in the predictions using the Spalart- Almaras and the Menter SST Turbulence 
Models 


¢ Additional effort is necessary to understand and resolve the CF prediction differences 





The conclusions of this study are shown in the Figure. 
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OVERFLOW. Facts on Friction 


Scott Lawrence 
NASA Ames Research Center 


Annual Airframe Technical Review 
Anaheim, California 
February 9-13, 1999 


NASA Ames Research Center 





This paper presents results of three minor studies into the behavior of the 
OVERFLOW with respect to the prediction of skin friction drag on wing 
bodies at cruise Mach number and wind tunnel Reynolds number. The studies 
include a preliminary assessment of the behavior of the two new 2-equation 
turbulence models introduced with the latest version of OVERFLOW (v. 1.8f), 
an investigation into potential improvements in the matrix dissipation scheme 
currently implemented in OVERFLOW, and an analysis of the observed 
sensitivity of the code’s skin friction predictions to grid stretching at solid 
surface boundaries. 
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Current Understanding 


Skin friction prediction variations are 
a primary source of uncertainty in 
the CFD prediction of absolute 
drag levels in supersonic cruise. 


NASA Ames Research Center 
HSA ATAYFebruary, 1999/2 





This schematic describes what is believed to be the current understanding of 
the sources of CFD drag prediction uncertainties, at least for the prediction of 
cruise Mach number wing/body configuration performance. 


For “reasonable” grids and algorithms, pressure drag appears to be relatively 
consistently predicted, whereas skin friction drag can be quite variable. The 
factors that might potentially cause this variability are shown as: code, grid, 
algorithm, and turbulence model. Of these factors, the code itself is not 
thought to cause uncertainty. Whereas poor choices of grid and algorithm have 
been found to produce a wide range of skin friction drag estimates, it is the 
author’s belief that the processes currently in use employ choices for these 
parameters that are suitable to reduce the uncertainty in friction drag 
predictions to acceptable levels (< 1 count) for wing/body configurations. For 
these parameters, poor choices can, to some extent, be identified without 
experimental comparison. For example, using artificial dissipation at the wall 
gives different results than scaling dissipation down near the wall; in the 
absence of experimental data, it 1s reasonable to conclude that the latter results 
are the more accurate. With a few exceptions (e.g., the low Reynolds number 
laminar flow produced by the Spalart-Allmaras model) the same cannot be 
said for the variation associated with the use of different turbulence models. 
Thus, the choice of turbulence model appears to have become the primary 
source of uncertainty in CFD drag predictions. 
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Outline 
e New OVERFLOW turbulence models 


e Grid-to-Grid comparisons 


¢ Matrix dissipation in OVERFLOW 


¢ Grid stretching analysis 


¢ Summary/Conclusions 


N ASA Ames Research Canter 
HSR ATR/February, 1999/3 





The outline of the presentation is shown here. The structure is intended to 
address elements of the previous figure, namely turbulence model, grid, and 
algorithm effects on friction drag. 
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OVERFLOW (1.8f) Turbulence Models 


Flat Plate Test Case: Mach 2.4 Wind Tunnel Conditions 
0.004 
e New models in OVERFLOW vt1.8f: ——— B-L 
~ k-w (Wilcox/Huang) ee ss (93) 
- SST (Menter) 


¢ Local skin friction groups into two 
distinct levels at moderate Re, 


— B-L 


¢ Average skin friction complicated aa ae (93) 


by low-Re, behavior --~- SST 
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The version of OVERFLOW exercised in this study is v1.8f. This version 
includes two new (to OVERFLOW) 2-eq turbulence models: the k-w model 
developed by Wilcox and the SST model (Menter), which is a modification of 
the Wilcox model. 


The results shown here, for a Mach 2.4 flat plate, indicate that the two new 
models produce significantly different local skin friction levels at moderate to 
high Re. In fact, the k-w model tends to track the local skin friction produced 
by the Spalart-Allmaras (S-A) model, whereas the SST model tracks the 
Baldwin-Lomax (B-L) skin friction. 


At low Re, there is a wide variation observed between the four models used 
and this variation has a substantial effect on the behavior of the average skin 
friction distributions. In particular, the average skin friction generated by the 
S-A model is severely affected by the laminar run region at low Re. This 
results in a much lower slope with Reynolds number than any of the other 
models. 
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OVERFLOW Turbulence Models: Integrated Drag 


TCA Wing/Body: Mach 2.4 Wind Tunnel Conditions 
Single-zone C-O grid (modified LB grid), truncated (3148”) 
New models add uncertainty: 
¢ Spread of 5 counts across model families 


° k-@ model notoriously sensitive to freestream turbulence 


Cor 


Counts 


SA92 SA93 SAQ3t 


NASA Ames Research Center 
HSA ATA/February, 1999/5 





This bar chart indicates results of the four turbulence models (plus two 
modifications of the S-A model) in application to the TCA wing/body. The 
grid used is a single-zone C-O topology grid truncated at the sting location 
(x=3148”). The dimensions of the grid are 97x241x65 and the wall spacing 
provides an average y+ at the surface of approximately 0.7. Stretching was 
applied beginning at the third point off the surface. 


The B-L and S-A (as implemented in OVERFLOW, SA93) models show 
approximately two counts of variation. This is caused by the presence of the 
laminar flow region near the wing leading edge over most of the span. 
Application of the original, or classic, S-A model (SA92), produces fully 
turbulent flow from the leading edge and similar integrated friction drag to that 
obtained with the B-L model. Overtripping, as suggested by Spalart, can also 
be used to cancel the effect of the laminar flow, though this process requires 
the specification of trip lines within the CFD grid. In addition, results are quite 
sensitive to a length parameter that must be specified by the user, which limits 
the usefullness of this approach as a predictive tool. 


Note that the variation with turbulence model is increased by virtually an 
order-of-magnitude with the addition of the 2-eq. models. It should be noted 
that the k-w model is known to be sensitive to freestream turbulence levels 
specified in the code. This effect has not been adequately investigated here. 
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OVERFLOW Turbulence Models: Local C, 


TCA Wing/Body: Mach 2.4 Wind Tunnel Canditions, Singie-Zone Grid 


NASA Ames Reasearch Canter 
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In an attempt to localize the sources of the friction drag differences, local skin 
friction has been plotted at four equally spaced span stations. Styles are 

applied in a manner similar to the flat plate case. It is apparent that the local 
skin friction on the lower surface (higher levels) behaves somewhat flat plate- 
like in a relative sense. That is, the S-A and k-w models produce similar 
distributions and the B-L and SST produce similar distributions, but lower than 
those of the other models. The behavior on the upper surface is more 
complicated, indicating a variation in the response of the various models to 
pressure gradient. Also, the outboard section (z=550’) is complicated by low 
Re effects. 
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OVERFLOW Grid-to-Grid: Integrated Drag 


TCA Wing/Bady: Mach 2.4 Wind Tunnel Conditicns 


FS ae ee eee Total Skin Friction Drag 


@ SingleZone 
.|@Overset 


B-L SAY2 SA93 SA93t SST 


Overset results due to Chaney, BCAG 


id 


SA93t - SA (93 mod) with “over-tripping 


END MADDEN RR 





HSA ATRFebruary, 19997 


Skin friction drag predictions for the TCA wing/body are shown here for 12 
different solution processes. Variables include turbulence model, code, and 
erid. For each turbulence model, variations tend to be modest (< | count of 
frag), except for the SST model which indicates a significant sensitivity to 
grid. Overset grid solutions were computed by Chaney at BCAG. 
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OVERFLOW Grid-to-Grid: SPAN = 100” 


TCA Wing/Body: Mach 2.4 Wind Tunnel Conditions, Singie-Zone Grid vs. Overset 
SA93 


¢ Cut through coliar/wing overlap 


¢ SA93 “tripped” by fuselage flow 


¢ Grid effects slightly greater for SST 


N ASA Ames Research Center 
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In order to localize the sensitivity to grid observed with the SST model, local 
skin friction coefficients are plotted at the four span stations previously used. 


Results are shown here at z=100” for three turbulence models, including the 
SST model. Some chatter is observed in the overset grid results because this 
Station passes through the wing/collar overlap region. Also of interest is the 
lack of laminar flow observed in the SA93 results. This is caused by the 
proximity of the fuselage boundary layer, which provides an influx of 
“turbulence” sufficient to numerically “trip” the wing boundary layer. Some 
grid sensitivity is observed in all three plots, but slightly more sensitivity is 
apparent in the SST distributions. 
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OVERFLOW Grid-to-Grid: SPAN = 250° 


TCA Wing/Body: Mach 2.4 Wind Tunnel Canditions, Single-Zane Grid vs. Overset 
SA93 


== Sl ngle-zone 
--~- Overset 


¢ Middle of inboard wing 


¢ SA93 shows laminar region 





¢ Significant grid effects only w/ SST 


| —— Single-zone Single-zone 
Overset ---- Overset 
> UPS 


4 
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At z=250”, the grid sensitivity of the SST model appears to be more severe 
than at the inboard station. The overset result produces high skin friction on 
both the upper and lower surfaces. Results for the two versions of the S-A 
model appear to be insensitive to grid and code at this span location. The SA93 
model exhibits the characteristic undershoot or laminar flow region near the 
leading edge on both upper and lower surfaces at this location. 
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OVERFLOW Grid-to-Grid: SPAN = 400” 


TCA Wing/Body: Mach 2.4 Wind Tunnel Conditions, Single-Zone Grid vs. Overset 


* Slightly inboard of wing break 


¢ Reduced grid effects on aft lower 
surface (SST model) 


* SA92 and SA93 show good grid- 
to-grid and code-to-code 
agreement 
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Near the outer edge of the inboard wing section, the grid sensitivity of the SST 
model is more modest, especially on the lower surface. Again, little grid or 
code sensitivity is observed with either S-A model. 
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OVERFLOW Grid-to-Grid: SPAN = 550° 
TCA Wing/Body: Mach 2.4 Wind Tunnel Conditions, Single-Zone Grid vs. Overset 
SA93 


¢ Outboard of wing break —— Single-zone 
-~-- Overset 


° Flow nearly 2-D ee 
¢ SST model appears nearly grid 
independent 








« Some code dependence showing 
w/ SA93 


Single-zone i \3 —— Single-zone 
---- Overset ; --—- Overset 
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Finally, on the outboard wing, the grid sensitivity of the SST model is nearly 
gone. Some sensitivity to code (or grid) is observed with the SA93 model, but 
it is of a relatively small magnitude and is limited to a relatively small area. 
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OVERFLOW Dissipation Dilemma 


(Faster, Better, AND Cheaper?) 


Mach 2.4 Fiat Plate 


Convergenc 
*| Profiles 


‘|Ct 





Single-Zone TCA 


SE ee rae 
/—— Sealar (SST); 


---- Roe (SA93) | 
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The OVERFLOW code provides essentially three algorithm dissipation 
options. In terms of three criteria of interest in the HSR program, none of the 
algorithm options is entirely satisfactory. Each has a critical drawback: the 
scalar dissipation method produces velocity profiles with overshoots near the 
boundary-layer edge, the upwind method has convergence problems in many 
cases, and the matrix dissipation tends to produce skin friction significantly 
lower than the other methods. 


The matrix dissipation scheme is similar to the scalar method in that a 
second/fourth dissipation operator is applied. However, it is applied to a 
difference of flux-like terms rather than conserved variable-like terms. With 
the scalar dissipation method, the user is given the option of modulating the 
dissipation near the wall to lessen its impact on the flow at the wall and, thus, 
the skin friction. This has been found to be critical to obtaining accurate 
friction drag estimates from OVERFLOW. In the present version of 
OVERFLOW, the matrix dissipation scheme has no such option. 
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MODIFIED Matrix Dissipation 
Mach 2.4 Flat Plate 


SMOO = 1 option added (analogous 
to Scalar method) 


- Upwind-like behavior for Cf and 
profiles 


0.004 


- Convergence behavior uncertain 





Hteration 
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If an option is implemented for the matrix dissipation method analogously to 
what is done for the scalar dissipation, results are obtained as shown here. Skin 
friction results are observed to align with the scalar and Roe results, and 
velocity profiles are obtained which are virtually indistinguishable from those 
obtained using Roe. Convergence histories with the modified matrix method 
(Matrix-1) are similar to those obtained with the original matrix scheme. 
However, the convergence histories for the flat plate problem don’t provide a 
good example of the problems that can be experienced with the Roe scheme. 


Preliminary calculations using Matrix-1 on the TCA wing/body indicate that 
the skin friction and profiles are improved in that case as well; however, it was 
observed to be somewhat less robust than the scalar method and the Roe 
method, so further work is required. 
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OVERFLOW Pseudo-Finite-Volume 


Effect of Valume Averaging on Grid Sensitivity 
simple Flat Plate Example (neglecting streamwise averagin 
For stress at secondary point 2: 
*) _ Uy ~My 
Oy}, Ay 


OVERFLOW: Ay* = (Ay,.+Ay,.)/2 = (Ay,+2Ay,t+Ay,)/4 
Ay* = Ay+(y-1) Ay/4 


For stress at secondary point 3: 
OVERFLOW: 


Ay* = yAy+(y-1)? Ay/4 


There is a sensitivity to stretching (y) everywhere, but eta. patne 
especially at points where there is a sharp change in jJesudescee antes 
stretching. 
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Finally, the sensitivity of OVERFLOW skin friction drag predictions to grid 
stretching at the wall has been known for several years. However, the 
underlying cause of the sensitivity has not been clearly understood. In attempt 
to provide a better understanding of the problem, a simple analysis was 
undertaken using a flat plate example. 


The pseudo-finite-volume approach taken by OVERFLOW applies flux 
conservation to implied volumes (shown as dashed red line) for which 
geometry information is determined using averages of the geometries of the 
surrounding “true” volumes. As a result, for the evaluation of the shear stress 
at the pseudo-cell face, an average of Ay from the surrounding grid cells is 
effectively used rather than the actual Ay that would be used in a more formal 
finite-volume approach. The error introduced through this approximation is 
proportional to (y-1)? in regions where the growth factor, y, is constant. This is 
a relatively modest error as long as y is small (1 <y< 1.2). 


The real problem occurs at the wall where the grid is effectively reflected to 
get boundary information. As such, y= | at the wall, which introduces a 
discontinuity in y, and the error in Ay incurred by averaging increases by a 
factor of 1/(y-1). The error can be reduced again by extrapolating into the wall 
rather than reflecting, but the error is not completely eliminated without 
changing the way the fluxes are evaluated. It is not clear at this time, to what 
extent the code needs to be modified to address this problem. 
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Conclusions 


OVERFLOW’s new turbulence models add to existing 
uncertainty in CFD friction drag predictions 
~ models produce significantly different flat-plate C, 
distributions 
experimental comparison (local or integrated) 
— two-equation model(s) appear much more sensitive to grid 
than Spalart-Allmaras (one-eq.) 
“controlled overset/single grid study 


Matrix dissipation can probably be modified (ala 
scalar) to produce “good” skin friction. 


Jacobian averaging is believed to be the source of 
OVERFLOW ’'s sensitivity to grid stretching. 


“move averaging into flux routines 
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Conclusions are drawn from the work presented. 


Uncertainty appears to be increased by the new turbulence models available in 
the new version of OVERFLOW. The uncertainty stems from differences in 
flat plate skin friction levels, as well as some observed grid sensitivities with 
the SST model. Presumably, the k-w model would exhibit similar tendencies if 
applied with the grids used here. It is felt that some of the differences in flat 
plate skin friction predictions could be sorted out through the use of 
experimental data, such as existing flat plate local skin friction and/or 
integrated skin friction through the symmetric model test. The source of the 
grid sensitivity of the 2-eq. models could be identified with controlled grid 
studies involving overset and single-zone grids on a simple geometry. 


It is believed that modifications to the matrix dissipation method currently 
implemented in OVERFLOW could be effectively applied to correct skin 
friction underpredictions, though further work is needed to establish exactly 
how the modifications could be applied to maximize the robustness of the 
resulting algorithm. 


Finally, simple analysis indicates that the sensitivity displayed by 
OVERFLOW to grid stretching at the wall is caused by the volume averaging 
used. The solution to the problem appears to involve moving the averaging 
into the flux routines to provide “flux-specific” averaging, though the 
performance penalty that would be incurred is unknown. 
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Canard Integration Wind-Tunnel Tests 
and Computational Results 


Todd E. Magee, Paul Kubiatko, Suk C. Kim 
The Boeing Company, Phantom Works 
Long Beach, California 90807 


Hamid Jafroudi 
Alpha STAR Corporation 
Long Beach, California 90807 


Abstract 


This paper presents experimental and CFD investigations into the effect of canard 
integration on the TCA (Technology Concept Aircraft) aerodynamic performance and 
Stability & Control. Specifically, results from the Supersonic Canard Integration Wind 
Tunnel Test (Test 1705) at the NASA Langley UPWT (Unitary Plan Wind Tunnel) and 
the Transonic Canard Integration Wind Tunnel Test (Test 508) at the NASA Langley 16- 
foot Transonic Wind Tunnel are presented. Comparisons of this test data with 
extensive CFL3D computational fluid dynamic (CFD) solutions are shown, as well as 
the results of a CFL3D directional stability & control study for the full TCA configuration 
(wing/body/canard/empennage). CFL3D shows good correlation with the test data for 
all configurations tested, except the ACC (alternate controls concept) canard 
configuration at a canard incidence angle of 4°. Navier-Stokes results for both the 
Baldwin-Lomax and Menter’s x-@ SST models are also shown. The CFL3D directional 
stability study indicates that the interaction of the PTC (Preliminary Technology 
Concept) canard tip vortices with the tail of the aircraft have a significant effect on 
directional stability. This tip vortex interaction is similar to low speed results which 
indicate a strong correlation between canard tip height and the angle-of-attack at which 
representative HSCT configurations become directionally unstable. Conditions for the 
studies presented in this paper range from M,=0.6, 0.9, and 1.1 for the transonic studies 
to M.=2.4 for the supersonic studies. Reynolds numbers for all studies were 
approximately 4 million/ft. 
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Run 115, Pt 5100, no canard 
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CFL3D Euler, M..= 0.9, « = 3.5° , B = 3.0° 
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low—mount canard (i = 4°) 


high—mount canard (i = 4°) 
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Unstructured Navier-Stokes Analysis of Wind-Tunnel Aeroelastic 
Effects on TCA Model 2 


A method is presented which accounts for aeroelastic effects in Navier-Stokes 
computations of low aspect-ratio wind-tunnel models. It is applied and 
validated in this presentation on the TCA Model 2 configuration. 
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Configuration Aerodynamics Technology Development 
Session 2: Analysis and Viscous Drag Prediction 


Goals Demonstrate Significant L/Dmax Gains 


Objectives Robust Analysis / Realistic Aerodynamic Efficient Engine / 
Testing Methods Design Optimization Airframe Integration 


a 
Challenges | Validation; | Viscous Effects F Multipoint Conditions | | Power Effects 


Analytic Methods ; Design Test Programs 
mppC anes and Applications Development and Techniques 


Methods Down Select Nace lle / Diverter 7 WT Database 
Design Integration 
I 


Program Viscous Drag Prediction WT Data Corrections 
; —— Tech. Baseline ; 
Cruise Point Optimization Development High Rel No. Testing 
Multi-Point Optimization Aero S&C PIE Test Program 
Development 


S&C CFD Predictions 
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e The Problem 

¢ Objective 

¢ Approach 

¢ Methodology & Procedures 
¢ Results 

e Summary 
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Problem - WT Aeroelastic Effects 


For rigid low AR HSR 
planforms, CFD predicts: 
¢ Lift slightly high 

¢ PM significantly low 


© Experiment 
—— Rigid CFD 
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The Problem: Wind-Tunnel Aeroelastic Effects 


The issue to be addressed is how to account for the effects of wind-tunnel 
model deformation in numerical computations of low-aspect ratio planforms 
without relying on measured model deformation data. 


During wind-tunnel tests of the low aspect ratio planforms characteristic of 
HSR configurations, the load induced displacement and twist of the thin outer 
wing panels tends to unload those panels. This results in a decrease in lift and 
increase in pitching moment. The variation of drag with lift is not adversely 
affected. Thus, it is important to include the model deformation within 
numerical computations to remove the noted aeroelastic-induced 
inconsistencies. 
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Objective 


¢ To improve CFD correlations with TCA Model 2 
UPWT data by coupling a simple calibrated structural 
model* into Navier-Stokes computations of wind- 
tunnel configuration 


¢ Validate over a range of experimental test points 
using a single point calibration of structural model 


* This approach has been successfully applied to 
high aspect-ratio Advanced Subsonic Transport 
configurations 





Objective 


The aim of this work is to demonstrate a simple technique which accounts for 
aeroelastic deformations experienced by HSR wind-tunnel models within CFD 
computations. With improved correlations, CFD can become a more effective 
tool for augmenting the post-test understanding of experimental data. 


The present technique involves the loose coupling of a low-level structural 
representation within the ELAPS code, to an unstructured Navier-Stokes flow 
solver, USM3Dns. The ELAPS model is initially calibrated against bending 
characteristics of the wind-tunnel model. 


The strength of this method is that, with a single point calibration of a simple 
structural representation, the static aeroelastic effects can be accounted for in 
CFD calculations across a range of test conditions. No prior knowledge of the 
mode] deformation during the wind-on test is required. 


This approach has been successfully applied to the high aspect-ratio planforms 
of subsonic transports. The current challenge is to adapt the procedure to low 
aspect-ratio planforms typical of HSR configurations. 
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Approach 


Set up 6-plate ELAPS structural model of TCA Model 2 
— Input geometric and material characteristics to ELAPS 
— Perform pre-test bending calibration of wind-tunnel mode! 


Compute USM3Dns/ELAPS coupled N-S solution of TCA 
— WI/B/N/DE at Mach 2.4, a=3.5°, Re,=6.4M 


— Manually move viscous grid with POSTGRID and restart sol’n 

— Integrate F&M for W/B/N/D 
Compare a=3.5° solution with Pitching Moment data 

— Scale calibrated wing deformation until PM data is matched 
Run additional AOA’s with final calibration and validate 
with other test point data 

— F&M and wing deflectionAwist measurements 





Approach 


This slide outlines the overall approach. Additional details will be described in 
the remaining slides. 


The approach is to construct a 6-plate structural representation of the TCA 
Model 2 with the Equivalent Laminated Plate Solution (ELAPS) code. The 
ELAPS model requires some geometric and material characteristics to be 
prescribed as input. This model is calibrated by a pre-test model bending 
procedure. 


A coupled USM3Dns/ELAPS Navier-Stokes solution will be generated on the 
full W/B/N/D/E configuration at the cruise point, Mach 2.4 and angle of attack 
3.5 degrees. Force and moments will be integrated only over the W/B/N/D 
portion for comparison with the Model 2 data. Grid movement will be 
performed manually using the POSTGRID code, which is part of the 
VGRIDns unstructured grid generation system. 


The o=3.5° solution will be compared with the experimental pitching moment 
data in order to fine tune the ELAPS calibration. Once the calibration is 
complete, the method will be tested at other angles of attack to validate the 
generality of the procedure. Correlations will be made against experimental 
force/moment and wing deflection/twist data. 
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The TetrUSS System 
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The TetrUSS System 


The flow computations are performed using the Tetrahedral Unstructured 
Software System (TetrUSS). This is a modular set of codes for solving Euler 
and Navier-Stokes problems on complex configurations. The system is based 
on tetrahedral volume elements which are extremely flexible for modeling the 
most complex of configurations. 


TetrUSS consists of a geometry setup tool, GridTool, tetrahedral grid 
generator, VGRIDns/POSTGRID, flow solver, USM3Dns, and graphical 
analysis tool, VIGPLOT. These tools been relatively easy to use by a broad 
range of users. Because of it’s loose coupling, other codes have have also been 
utilized within the framework by various users as denoted by parenthesis. 


The TetrUSS system supports additional engineering capabilities such as 
interacting boundary layer (Euler/IBL), aeroelasticity, iterative design, and 
propulsion effects. These capabilities will be described in more detail in the 
next slide. 
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Modular TetrUSS Capabilities 


* Can substitute other flow solvers 





Modular TetrUSS Capabilities 


TetrUSS includes a modular capability for computing aeroelastic effects, 
iterative design, and/or interactive boundary layer. Each component is 
maintained independently and is coupled to the system by pre- and post- 
processing utilities. The capabilities can all be used simultaneously, or in any 
combination. 


During the script-driven update cycle, a solution file is written by the flow 
solver and is used by the sub-modules. Both the aeroelastic and design features 
involve grid movement, whereas the IBL feature generates transpiration 
boundary condition velocities. 


The aeroelastic component converts the loads from the flow file into input for 
the structural code. The present options for structural modeling include a 
coupling with the ENSAERO system and the ELAPS code. Higher-level Finite 
Element Model (FEM) structural representations can be input through 
ENSAERO via modal or stiffness matrices. 


While higher-order structural modeling is an option of this system, the present 
work is focused on utilizing a low-order structural representation through a 
simple plate model in the ELAPS code. 
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USM3Dns - Salient Features 


- Tetrahedral cell-centered, finite volume 


Euler and Navier-Stokes 
— Spalart-Allmaras turbulence model with wall function 


Time Integration - Implicit GS and Explicit RK 


Roe’s upwind FDS with flux limiting 


Platforms 
~ Cray vector machines with multitasking 
-— Massively parallel computers or loosely coupled 
workstations 
Coming capabilities 
— Advanced 2-eqn turbulence models (in testing) 
— Low-Mach preconditioning (in testing) 
— Time-accurate moving grid (planned) 





USM3Dns - Salient Features 


USM3Dns is a tetrahedral cell-centered, finite volume Euler and Navier- 
Stokes flow solver. Turbulence is modeled by the Spalart-Allmaras one- 
equation model. The sublayer region of the boundary layer can also be 
modeled analytically by a wall function in order to reduce near-wall grid 
density. 

Solutions are advanced in time by either an implicit Gauss-Seidel scheme or 
the Jameson Runge-Kutta explicit time-stepping. Cell flux is computed with 
the Roe’s Flux Difference Splitting (FDS), which can be limited by either a 
MinMod or SuperBee limiters. 


USM3Dns will run on Cray vector processors with multitasking efficiencies of 
6 out of 10 processors. Another Message Passing Interface (MPI) version of 
the code will run on massively parallel computers or loosely coupled 
workstations. Superlinear performance has been demonstrated on an Origin 
2000 with 150 processors. 

Additional capabilities such as 2-equation k-epsilon and Reynolds stress 
turbulence models, and low-Mach preconditioning have been implemented and 
are currently being tested in a non-release version of USM3Dns. The current 
plan is to a add time-accurate moving grid capability. 
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ELAPS Structural Code 


e ELAPS - Equivalent Laminated Plate Solution 
-— Computationally efficient method based on plate model 
— Employs Rayleigh-Ritz method with a global displacement 
polynomial function 
e Less cumbersome to use than FEM 
— Suitable for preliminary aeroelastic effects 
e Easily coupled with USM3D 
— Deformed surface points displaced vertically 
— Volume grid moved manually* (at this time) via. POSTGRID 


* Can be easily automated in the future 





ELAPS Structural Code 


The ELAPS code, developed by Dr. Gary Giles at NASA LaRC, is based on an 
equivalent laminated plate technology. For the present work, it can be applied 
to set up a computationally efficient polynomial plate representation of the 
wing planform. The Rayleigh-Ritz method is used to generate a polynomial 
function of vertical displacement over the wing planform. 


While one could utilize a more complex finite-element model in the TetrUSS 
system, the advantage of the present approach is in its simplicity and rapid set 
up time. 

Grid movement involves moving both the thin-layered ‘viscous’ tetrahedra as 
well as the outer ‘inviscid’ cells. This function is performed manually with the 
POSTGRID code at this time, but could be easily automated in the future. 
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Twist Calibration Procedure for Structural 
Plate Representation of Wind-Tunnel Model 





Twist Calibration Procedure 


This slide describes a simple pre-test calibration procedure of the ELAPS 
structural representation. Weights are hung a selected points on the outer wing 
panel of the wind-tunnel model. Measurements of model deformation are 
recorded at several locations on the configuration. 


A simple solid plate representation is constructed of the model structure within 
the ELAPS code. Preliminary estimates of material properties are also 
prescribed. Simulated point loads are applied to the ELAPS model which 
produce an initial “computed” twist distribution which is compared to the 
measured values. The material properties are adjusted within ELAPS until 
there is reasonable agreement between the twist deformations from the 
computed and measured data. 


Once this calibration is achieved, the coupled ELAPS model is then included 
in the Navier-Stokes computations across the test polar. 
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Point-Load Calibration of TCA Model 2b in UPWT 


© Ft. Load Cal. 
—— Calibrated ELAPS 
---—-- Uncalibrated ELAPS 


Load Point O\_, 
101.46 1b 


2.0 
0.0 0.2 04 06 0.8 1.00.0 0.2 04 06 08 1.0 
2y/b 2y/b 





Point-Load Calibration of TCA Model 2b in UPWT 


This slide illustrates how the weights are applied to the inverted TCA Model 
2b in the Unitary Plan Wind Tunnel. The inset photo depicts the recording of a 
measurement of tip deflection. 


The plots contain the measured twist distribution induced by point loads 
applied at the 72-percent and 87-percent span stations. The dashed line is the 
raw “uncalibrated” ELAPS twist distribution based on the prescribed 
geometric and material properties of the model. The solid “calibrated” line 
was derived by adjusting the material properties in the 6-plate ELAPS model 
until a reasonable approximation of the measured data was obtained. As can be 
observed in the plots, some discrepancies do remain particularly at the tip with 
the inboard load point. 
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Changes to ELAPS Calibration Strategy 


¢ Initial N-S calculations with standard pre-test calibration 
procedure achieved about a 40-percent improvement in 
pitching moment correlation at a=3.5° 


¢ Adjusted calibration of 6-plate ELAPS model by 


doubling its surface displacements 
— achieved close match with pitching moment data at a=3.5° 


¢ Found last minute oversight in calibration 
— Bending calibration done on Model 2b from T1846 (solid wing) 
— Correlations made with Model 2a from T1671 (has flap cutouts) 
— 11671 measured deformations twice that of T1846 





Changes to the ELAPS Calibration Strategy 


After the initial Navier-Stokes calculation using the calibrated ELAPS model, 
the correlation with the experimental pitching moment had improved only 40- 
percent. It was initially assumed that the robustness of the 6-plate ELAPS 
model was not adequate to yield a straight forward calibration for low-aspect 
ratio planforms. A further adjustment was then made by doubling the surface 
displacements produced by the ELAPS model. This yielded a close match with 
the experimental pitching moment data at a@=3.5°. 


The source of this shortcoming was later determined to be an oversight on 
which Model 2 was calibrated. The bending calibration was performed on the 
Mode] 2b from T1846 which had a solid wing structure, whereas the 
experimental data was taken from Model 2a of T1671 which was considerably 
more flexible due to leading- and trailing-edge flap cutouts. After a last minute 
check of optically measured model deformation data from the two tests, it was 
determined that twist measurements differed normally by a factor of two. 
Additional calculations will be made after the workshop to resolve this 
inconsistency. 
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Solution Characteristics 


e Thin-Layer Tetrahedral Grid 
— 1,428,427 cells (897,078 in ‘viscous’ layer) 
— Average y* of 1st node is 20.1, and cell centroid is 4.5 
— Approx. 18 tetrahedra across mid-chord BL 


e Flow conditions | 
—- M=2.4, o=3,3.5,4,&5deg, Re,=6.4x108 


e Flow model 
— Full Reynolds-averaged Navier-Stokes 
— Spalart-Allmaras turbulence model 
— Wall function below first node (y* approx. 20) 





Solution Characteristics 


A tetrahedral grid of 1.4 million cells was generated with VGRIDns. The 
average y* obtained by averaging the values at the first node above the surface 
on all viscous surfaces was 20.1. Noting that the distance of the centroid of a 
first layer tetrahedral viscous cell above the surface is 1/4 that of the first node, 
the average y* at the cell centroids was 4.5. 

The flow conditions are Mach 2.4, a=3, 3.5, 4, & 5 degrees, Re,=6.4 million. 
The computations incorporate the full Reynolds averaged Navier-Stokes 
equations above the first layer of tetrahedral cells using the Spalart-Allmaras 
turbulence model, and model the flow within the first layer (below y* approx. 
20) using a wall function. 
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HSR Technology Concept Airplane 
Tetrahedral Viscous Grid (1,428,427 cells) 
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HSR Technology Concept Airplane 


This slide depicts the surface and symmetry plane grid used in the 
computation. Note the clustering of the ‘viscous’ cells along the 
fuselage/symmetry plane boundary. All solutions were computed on the full 
W/B/N/D/E configuration, but forces and moments were integrated only over 
the truncated W/B/N/D geometry for comparison with supersonic wind-tunnel 
data. 

Regarding movement of the nacelle/diverters under aeroelastic deformation, 
no special treatment was applied. TetrUSS has a capability for slaving the 
movement of components. But this could not utilized for TCA due to the tight 
integration of nacelle/diverter with the wing. Fortunately, the required nacelle 
movement was very small. 
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Solution Strategy 


Run 700 cycles on initial “rigid” grid 

— Run additional 500 cycles to generate “rigid” reference sol'n 
Move grid using preceding flow solution 

~— Move surface grid front with PREELAPS utility 

— Move volume grid with POSTGRID 

Restart USM3Dns from rigid solution for another 100 
to 150 cycles with USM3Dns 


— repeat process 3 times to converge aeroelastic solution 


Resource requirements 
— Memory: 250MW 


— CPU time: 12.5 C90 hours per AOA for initial 700 cycles; 
7 hours to converge aeroelastic solution 





Solution Strategy 


The aeroelastic solutions were generated through a series of restart runs 
beginning with a partially converged ‘rigid’ solution, and continuing after each 
movement of the grid. The grid was moved manually with POSTGRID. A total 
of three movements/restarts were required to converge the aeroelastic 
deformation at each angle of attack. Each case required 250 megawords of 
memory and 20 hours of Cray C-90 time. 
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Effect of Calibrated USM3Dns/ELAPS 
on F&M Correlation 


© Experiment 
--}--- USM3Dns: Rigid 
TCA Model 2 --w-~ USM3Dns: ELAPS 
M=2.4, Re,=6.4M wnee- COPLED (Deformed) 


a Calibration Point 
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Effect of Calibrated USM3Dns/ELAPS on F&M Correlation 


This slide presents the CFD correlations with force and moment data from 
UPWT T1671. Note that the ‘rigid’ USM3Dns result over-predicts lift and 
significantly under-predicts pitching moment. The drag coefficient is in good 
agreement and includes corrections for nacelle base pressure and for trip drag 
by the k? method. 


The ELAPS calibration point is denoted by the solid triangle on the pitching- 
moment curve. The remaining open triangles were computed with USM3Dns 
using the calibrated ELAPS model. Note the good agreement of both the lift 
curve and the pitching moment with experimental data. Note further the 
excellent agreement with the CFL3D computation presented by Kuruvila, et al. 
The CFL3D result was computed on a deformed structured grid which 
matched the optically measured shape of the wind tunnel model. 


Thus, the present approach offers a simplified technique for addressing 
aerodynamic deformation effects of low-aspect ratio wind-tunnel models in 
CFD computations without prior knowledge of those deformations. 


References 


Kuruvila, G., Hartwich, P. M., Baker, M. L.: “Effect of Aeroelasticity on the 
Aerodynamic Performance of the TCA”, 1998 NASA High-Speed Research 
Program Aerodynamic Performance Workshop, NASA CDCP-1006, Volume 
1, Part 2, pp. 1589-1648 
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Twist Deformation on TCA Model 2 W/B/N/D 
UTWT Run 1671/85, M=2.4, Re.=6.4x106 


© Data 
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Twist Deformation on TCA Model 2 W/B/N/D 


This slide conveys the correlation of computed and experimentally measured 
changes to spanwise distribution of wing twist at the four angles of attack. The 
general agreement is good, but the computed curves do not reflect the 
inflection present in the experimental curves. In particular, the changes in twist 
are under-predicted in the mid-span region and over-predicted near the tip. 
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Vertical Wing Displacement on TCA Model 2 W/B/N/D 
UTWT Run 1671/85, M=2.4, Re=6.4x106 
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Vertical Wing Displacement on TCA Model 2 W/B/N/D 


A comparison of the experimentally measured and computed spanwise 
distribution of vertical displacements are presented here for the four angles of 
attack. The experimental displacements were measured optically at several 
spanwise locations on the wing surface. Time did not permit us to extract the 
computed displacements at the same locations, thus, the leading- and trailing- 
edge values are included on this figure. Note that the computed displacement 
curves tend to bracket the experimental ones, with the exception that the tip is 
over predicted. | 


The results on this and the prior two slides suggest that while the simple 6- 
plate structural representation does not resolve all of the details of the 
aeroelastic phenomenon, it does result in a good estimate of the aggregate 
aerodynamic performance loads on the configuration. 


A better detailed modeling might be achieved either by gaining more 
experience with the ELAPS model, or by applying a similar calibration 
procedure to a higher-order structural model. 
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Summary 


Applied simple point-calibrated structural model to 
include WT aeroelastic effects in N-S computations of 
low aspect-ratio HSR planform 


Demonstrated coupled system of unstructured N-S 
flow solver (USM3Dns) and low-level structural code 
(ELAPS) with movement of ‘viscous’ grid 


Improved N-S correlations with TCA Model 2 over 
range of test points 


— Force and moment test data 
— CFL3D solution with measured deformations 


Reasonable agreement of twist and displacement with 
data 





Summary 


The objective of this work was to demonstrate a simple technique which 
accounts for aeroelastic deformations experienced by HSR wind-tunnel 
models within CFD computations with a view toward improving the 
effectiveness of CFD as tool for augmenting the post-test understanding of 
experimental data. 


Such a method has been validated herein using an unstructured Navier-Stokes 
flow solver, USM3Dns, and a low-level polynomial plate rendition of the 
structure through the ELAPS code. Technical challenges such as modular 
coupling of codes and movement of ‘viscous’ grids on a complex geometries 
have been overcome. The supporting results demonstrate that the aggregate 
aerodynamic performance loads are well predicted. While there is some 
deficiency in distributed quantities, the agreement of computed twist and 
displacement of the wing are in reasonable agreement with the data. 


In summary, the present approach offers a simplified technique for addressing 
aerodynamic deformation effects of low-aspect ratio wind-tunnel models in 
CFD computations without prior knowledge of those deformations. 
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AIRPLANE Nacelle/Diverter Grids and Solutions 
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Update Cos #2 
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SYN107-MB Warped Nacelle Design Comparison, Wing/Body/Nacelle/Diverter 


M=2.4, no internal or base nacelle forces, entire fuselage 


—O—_ Baseline TCA SCALAR ——>—— Ames 1-03 CUSP 

—<——__ Ames 1-03 SCALAR ——@®— Ames 6-08 CUSP 
0.11 
5 0.10 
0.09 





60.0 65.0 70.0 75.0 80.0 


CD (pressure) mM 
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SYN107-MB Warped Nacelle Design Comparison, Wing/Body 


=2.4, entire fuselage 


mes 1-03 CUSP 
mes 6-08 CUSP 
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S YN107-MB Warped Nacelle Design Comparison, Nacelle/Diverter 
M=2.4 


——¥——__ Ames 1-03 CUSP (no internal nac forces) 
——®—-__ Ames 6-08 CUSP (no internal nac forces) 





-10.0 
0. > 10. 15.0 20.0 


CD (pressure) 
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Comparison of the 6-08 Configuration with Deformed and Original Nacelles 
M=2.4, no internal or base nacelle forces, entire fuselage 


——>——— SYNI107-MB 1-03 SCALAR 
——@—— SYN107-MB 6-08 Deformed nacelles, SCALAR 
——S— __ SYN107-MB 6-08 Original nacelles, SCALAR 
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SYN107-MB Rigid Nacelle Design Comparison, Wing/Body/Nacelle/Diverter 


M=2.4, no internal or base nacelle forces, entire fuselage 


— > Ames 1-03, SCALAR 
--~-@®-- Ames 12-10 design, SCALAR 
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SYN107-MB Rigid Nacelle Design, Wing/Body Drag Comparison 
M=2.4, entire fuselage 


—*——_-__ Ames 1-03, SCALAR 
---~@-- Ames 12-10 design, SCALAR 
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SYN107-MB Rigid Nacelle Design, Nacelle/Diverter Drag Comparison 


=2.4, no internal nacelle forces 


———— Ames 1-03, SCALAR 
---@-- Ames 12-10 design, SCALAR 
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Figure 5: C, distribution at near wing tip station. Navier-Stokes calculations, M = 0.80, Re = 2.5 million 


, Aeroelastic solution 
—- — ~, Sohd geometry solution 
+ + +, Wind tunnel data 
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Figure 7: Cp distribution at near wing tip station. Navier-Stokes calculations, M = 0.80, Re = 2.4 million 
, Aeroelastic solution 

— - -, Solid geometry solution 

+ + +, Wind tunnel data 
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8a: span station z = 0.194 8b: span station z = 0.387 





8c: span station z = 0.581 . 8d: span station z = 0.775 


Figure 8: Typical Business Jet Configuration. Drag Minimization at Fixed Cy. 

Rigid Design, M = 0.82, Cr = 0.35, 133 Hicks-Henne variables. Spar Constraints Active. 
Rigid Analysis at Fixed Cr. 

~ ~— -, Initial Pressures 

————, Pressures After 20 Design Cycles. 
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Rigid Wing Design 
-~--+--- Baseline Wing Design 


«10° 


Stress 


0.0 0.2 0.4 0.6 0.8 


Span position 


Figure 10: Spanwise Stress Distribution for the Rear Spar. 
Comparison of the Rigid Design and the Baseline Design. 
Wing Alone Configuration. 

Rigid Design. Drag Minimization at Fixed Cy. 

Aeroelastic Analysis at Fixed Cy. 
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Rigid Design, Aeroelastic Analysis 
iam tra Baseline Design, Aeroelastic Analysis 


Load 





0.0 50.0 100.0 150.0 200.0 250.0 


Span position 


Figure 9: Spanwise Load Distribution. 

Comparison of the Rigid Design and the Baseline Design. 
Wing Alone Configuration. 

Rigid Design, Drag Minimization at Fixed Cy. 
Aeroelastic Analysis at Fixed Cy. 
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14a: span station z = 0.194 14b: span station z = 0.387 
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14c: span station z = 0.581 


14d: span station z = 0.775 


Figure 14: Typical Business Jet Wing Configuration. Drag + Stress Minimization at Fixed Cr. 
Aeroelastic Design with Stress Penalty Function. M = 0.82, Cr = 0.35 

133 Hicks-Henne variables. Spar Constraints Inactive. 

Aeroelastic Analysis at Fixed Cy. 

~ — ~, Initial Pressures 


, Pressures After 13 Design Cycles. 
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Aeroelastic Design, Aeroelastic Analysis 
echeaiia arias. Baseline Design, Aeroelastic Analysis 


Load 


0.0 50.0 100.0 150.0 200.0 


Span position 


Figure 15: Spanwise Load Distribution. 

Comparison of the Aeroelastic Design and the Rigid Design. 
Wing Alone Configuration. 

Aeroelstic Design, Drag + Stress Minimization at Fixed Cy. 
Aeroelastic Analysis at Fixed Cz. 
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Aeroelastic Wing Design 
Se aS Baseline Wing Design 


#10 


Stress 


0.0 0.2 0.4 0.6 


Span position 


Figure 16: Spanwise Stress Distribution for the Rear Spar. 
Comparison of the Aeroelastic Design and the Rigid Design. 


Wing Alone Configuration. 
Aeroelstic Design, Drag + Stress Minimization at Fixed Cy. 


Aeroelastic Analysis at Fixed Cy. 
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